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文章 编号 :1005-3085(2011)03-0365-10 


不 可 压 渗流 驱动 问题 混合 有 限 元 /间断 有 限 元 
耦合 格式 的 误差 分 析 * 


杨 继 明 
(湖南 工程 学 院 理学 院 ， 湘 潭 411104) 


摘 要 : 多 孔 介质 中 渗流 驱动 问题 与 环境 污染 和 油 藏 开采 等 问题 密切 相关 ， 是 当今 的 研究 热点 ， 对 具有 分 
子 扩散 和 弥散 效应 的 不 可 压 渗 流 驱 动 问题 ， 本 文 用 混合 有 限 元 / 间 疡 有 限 元 耦合 格式 来 求解 ， 即 
用 混合 有 限 元 方法 求解 压力 方程 ， 用 对 称 内 罚 间断 有 限 元 方法 逼近 浓度 方程 . 运用 比 剪 切 算 子 更 
为 便捷 的 归纳 假设 和 插值 投影 ， 导 出 了 先 验 hp 误差 估计 . 
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分 类 号 : AMS(2000) 65M12; 65M60 中 图 分 类 号 : 0241.82 文献 标识 码 : A 


1 引言 


我 们 考虑 下 列 不 可 压缩 的 多 孔 介质 渗流 驱动 问题 


V-u=—V- (a(c)Vp) =4q, (x,t) ENx J, 

o% +u-Ve—V-(D(u)Vc) = (€-c)q, (z,t)EQx J, 
(x,t) € Nx J, (1) 
( 


x,t) € 02x J, 


u-n=O0, 
D(u)Vc-n=0, 
c(zx, 0) = co(Z), TER, 
其 中 9 是 二 维 的 有 界 多 边 形 区 域 ，68 是 其 边界 , J = (0,T], n 表示 9Q 的 单位 外 法 向 量 ; u(x,t) 
为 流体 的 Darcy 速度 ，p(z,t) 为 流体 的 压力 ，c(x,t) 为 注入 流体 的 浓度 ，$(z) > 0 为 有 效 介 质 
LRE, Du) 表示 扩散 /弥散 张 量 . 


D(u) = dml + ful (cu (u) + œ (I — E(u))), 


RF Eu) ARF uA MLke, HG, j)eBA 


UjU; 


(E(u); , = que’ 


dm > 0 是 分 子 扩散 系数 ; a: > 0 和 a > 0 分 别 为 纵向 和 横向 的 弥散 系数 . 4 为 井 内 的 流 
量 ， q= qt 十 9 ， 其 中 qt = max(q, 0), q = min(q, 0). 假定 4 和 党 都 是 有 界 的 ， € 表 示 注 入 
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井 处 注入 的 浓度 (q > 0) 和 产 出 井 处 剩余 流体 的 浓度 (q < 0)， 系 数 a(c) > 0 有 界 ， aele) 是 一 致 
有 界 且 关于 c 是 Lipschitz 连续 的 . 

混合 有 限 元 方法 能 够 使 压力 和 Darcy 速 度 取 得 相同 的 收敛 阶 ， 被 广泛 应 用 于 多 孔 介 质问 题 
的 数值 模拟 [3. 间断 有 限 元 方法 是 一 类 非 协 调 元 BA， 因 为 它 有 着 很 好 的 物理 和 数值 特性 ， 在 
数值 模拟 中 很 受 欢 迎 . 文献 [5,6] 对 于 只 有 分 子 扩散 效应 的 可 压 渗流 驱动 问题 ， 采 用 了 混合 有 限 
元 /间断 有 限 元 看 合格 式 ， 在 分 析 中 ， 文 献 5] 引入 了 前 切 算 子 ， 而 文献 [6] 运用 了 归纳 假设 口 . 
不 过 他 们 的 分 析 仅 只 限于 没有 弥散 的 特殊 情况 ， 实 际 上 分 子 扩散 和 弥散 的 影响 都 是 存在 的 ， 
这 给 分 析 带 来 了 一 定 的 困难 . 文献 [7] 中， 对 同时 具有 分 子 扩散 和 弥散 效应 的 不 可 压 渗流 驱动 
问题 ， 利 用 前 切 算 子 导出 了 半 离 散 的 混合 有 限 元 /间断 有 限 元 方法 的 耦合 格式 的 误差 估计 . 本 
文 研究 了 同时 具有 分 子 扩散 和 弥散 效应 的 不 可 压 渗 流 驱 动 问 题 混合 有 限 元 /间断 有 限 元 耦合 格 
式 ， 并 采用 归纳 假设 进行 误差 分 析 . 这 样 避免 了 象 文 献 [5,7] 那样 选择 剪 切 算 子 中 正常 数 的 麻 
烦 ， 而 且 归 纳 假设 可 以 直接 由 离散 格式 和 已 得 误差 估计 式 证 得 . 





2 混合 有 限 元 /间断 有 限 元 耦合 格式 


设 士 是 Q 的 拟 一 致 网 格 剖 分 (三 角形 或 四 边 形 ). Th 为 所 有 内 部 边 的 集合 ，h = max hs 为 
最 大 单元 尺寸 . 对 s > 0， 定 义 下 列 间 断 Sobolev 空间 


H° (Th) = {ve L?( (Q): v|s € H (E), EeT}. 


WE; € Th, Ej € Th, y = OE, NOE; E Tr (HINER n EI EL; HSM). 定义 we HT), s > 
1/2 的 平均 和 跳跃 为 


1 
{v} = a(( le) ly + (le,)ly)> = (vle) ly- (vle) ly- 
间 汤 有 限 元 空间 定义 为 
= {ve LQ) :vlp E€ P(E), EE Th}, 
HH P, (E) 表示 在 如 上 总 的 度 小 于 等 于 7r 的 多 项 式 空间 . 定义 空间 
V = H(div;Q) = {u € (L7(Q))?, divu € L?(Q)}, 
V° = Ho(div; Q) = {u € H(div;Q), u-nlaq = 0}, 
W =L (Q). 
SV x WAT HD Tr EE T T E Ve(Tn) x We(Th) X k (k > 0) 阶 Raviart-Thomas 空 
间 (RT). 相应 地 V? 的 子 空间 为 VE (Tr) = Vela) Nn Ve. 
本 文采 用 通常 的 Sobolev AFR (-, -) AR ||- mo 特别 地 ，|| :lloo WF || I. ATR, 
在 积分 式 do 和/ .dt 中 省 去 了 dz 和 dt， 即 分 别 用 /,.(g = B,7,0) 和 帮 : 表 示 关 于 空间 变量 


的 积分 和 时 间 变 量 积分 ，K, Ki (i € N) 表示 一 类 与 六 了 和 大 无 关 但 可 能 与 问题 的 解 相关 的 正常 
数 . 这 些 常 数 在 不 同 的 场合 下 到 不 同 的 值 . 用 = 表示 一 个 可 以 任意 小 的 正常 数 . 
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对 任意 的 % e D (Th) ENIRE A B(w; c, Y) 和 线性 泛 函 Lle y): 


B(u; c, y) = 5 [ Dave: vy- D /pmvemm 


EcT, yen 
z Z pow.: .92}[a + D fe Vey + Jo (c, 4), 
TEIR EETh 


Ee) = [ (@- cay, 
其 中 j 表示 7 的 尺寸 ， 





Jo (c, is a c 
saw) = 5 = fiw 
AWM. o > 0 为 罚 参 数 ， 它 在 边 7Y 上 取 值 为 cv， Hee en 

对 压力 方程 ， 我 们 采用 混合 有 限 元 方法 求解 .对 浓度 方程 ， 我 们 采用 对 称 内 罚 间断 有 限 元 
方法 有 逼近 .那么 问题 (1) 的 混合 有 限 元 /间断 有 限 元 耦合 格式 为 ， 寻找 UE LO(I;VE(T)), Pe 
L°(J;We(Tn)), C € L*(J;D,(Th))» Efa 


(V-U,w)=(q w), Yw E Welt), (2) 
(a(CWU,») -(V-¥,P)=0, Yv € VÈT), (3) 
(05.0) +BU;C, 9) = L(C,4), YY E DAT), (4) 


HE a(c) = 1/a(c). 


3 ”插值 投影 和 归纳 假设 
设 钼 入 分 别 为 w 和 p 的 投影 满足 
(V-u,w)=(g,w), Yw EWT), (6) 
(a(C ù, v) —(V-v,p) =0, Yv EVP). (7) 


Wp=u—-u,o=u-U, n=p-p, r=p— P. 假定 c(.,0) = 0， 由 文献 [2,8] 可 以 得 到 如 
下 投影 误差 估计 
in(k+1,wE—1) gee lB!) 


lol + al sK X "elles, |2] <x = “2a Mle. 
EE 


其 中 为 RT 空间 的 阶 数 ， we 为 p 的 正则 阶 . 
AS HE SCHR [9,10], SE € Thn, w € HTa), ESEM ARES c, rhe ERAK 
WK, UR Ü EP, (E), SER RERE 有 


| s>0, (8) 





ake 6-1/2 
lw 一 Wls, OE < KS lulle, E, 8> z =+ ô, ô = 0, P (9) 
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WEA chi, c— CM O(c 一 避 /8 满 足 上 述 加 逼近 性 质 . S= eE, E= TCO, 并 


HL E(-,0) =0, We-C=E+C. 
在 分 析 过 程 中 ， 需 要 用 到 下 列 归纳 假设 (将 在 第 5 节 进 行 证 明 ) 


llollLw;L» 0) < L, 


E 


2 








Z2(Ji r=) 全 


其 中 工 为 某 个 正常 数 . 
4 RED 


MASK (6)-(7) 分 别 减 去 式 (2)-(3)， 得 


(V-o,w)=0, Vwe Wr.(T), 


(a(C)o,v) — (V-v,7) = ((a(C) — a(c))u, v), Vu Ee V(t). 


Rw=7,v=a, H 

(a(C)e, o) — ((a(C) — a(c))t, e) = 0. 
应 用 Halder 不 等 式 推出 

|(a(C) — a(O))%, o| < KIE + él - lll. 
注意 到 a(C) 是 有 界 的 ， 于 是 ， 得 到 了 压力 方程 的 误差 估计 


llall? < K (1&1? + ell). 
设 (p u, c) 是 问题 (1) 的 解 ， 它 满足 


(%4) T B(u;c, 4) = L(e, 4), vye D-(Tn), E 


注意 到 四 


e-o-e-o-1| -(C+ 日 ， 如 果 q>0, 
0, 如 果 q<0. 


由 于 
u- Ue rete 


用 式 (15) 减 去 式 (4)， 并 取消 = So, TA 


CE DD Aston bats 


(10) 


(11) 


(12) 


(13) 


(14) 


(15) 


=~ (4%, %2) + an (u) - DU))Ve-n} |% on D(u) ~ DV) Ve: v 


+O (Pw -POY nd- ee vey 


yen 


Djo vea? 
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sS og ME r(e S) 





EET, 
十 DVS. C+E cat 
2 n y ] 一 ford q +5 & 
= > VWED,(T,), ted. (16) 


MIT; (i = 1,2,.… 12) 进行 估计 ， 我 们 需要 用 到 下 列 不 等 式 . 注意 到 ， 当 及 足够 小 时 


1Zlleo < Nilloo,n + llolloo,n < M, 


|< < 
ot L2(J;L°(Q)) 一 


其 中 用 到 了 归纳 假设 (10)-(11)、|| 名 oo,n A || H luro) 的 有 界 性 加. 
因为 [c] = 0， 所 以 Ts = 0. 利用 Cauchy-Schwartz 不 等 式 ， 可 得 




















Je [a 
L2(J;L°(2)) — I Ot lL2(7;L~ (0)) Ot 


ri se + x. rm so) + x(a to? 


iral < eÉ} + avave), al se] EI + Kl, Tas xe. 
因为 
PO < dm + max(a,a)U} <M, |S) < lvl SI. < 


VERE E(x, 0) = 0， 对 了 关于 t 积 分 并 用 分 部 积分 法 ， 推 得 
LHE fo U)V- VE)( oll E poveva 
tf a EN 
eg [mot g [so 
xf > f (ver + KP +% 
五 ETh 


对 于 Ts 和 To， 利 用 分 部 积分 、Cauchy-Schwartz REX. oF 的 有 界 性 以 及 文献 [6] 中 的 
迹 不 等 式 引 理 3.1 和 道 不 等 式 引 理 3.2， 可 知 


| [ T+ To] < | E Pove: n} [2 | + lei x [Paws nhic+ al 





< KD rh ||VCC,t) ml? + eg (Et) EC t) +K f y7 rhv n| 
YETh 0 yern 
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+K f Jg Go+K [ DT 





= ‘ -1 2 
rth VC nl + Ka ( mja 5)" SO a) [vec] 











Elcomref D frer f E ELSE h) 
对 2 关于 上 积分 并 用 分 部 积分 ， 有 
t ro, r ray ð 
allg T atenono- 5 > T f [alk 
< IE (EC, t), EC, DIKE i 7 Z Ikeo? xf IEE) 
+e fF (Wen + UES). 
+r =Å S0 (u) - DU))Ve-n}{é] - 二 [{e DW) YG) 
yer,” 7 TEL 
-5 [ {Fw -Dw)ve-n}ie + ee (u) - DU))Ve- ve) 
yer, “7 
-PW ~ DUVE: Ve) + (D -DOVE ve) = Dh 








[ev < tol + lol, [EEP] < |] + |, 


根据 D(u) 的 定义 、c 和 Ve 的 有 界 性 、 带 6 的 Cauchy-Schwartz 不 等 式 以 及 文献 [6] 中 的 迹 不 等 
式 引 理 3.1 和 逆 不 等 式 引 理 3.2， 可 以 得 到 


[ nisk Y WG (FA) e-news + ST [ec]? 
< Koz? (o6 DI + lol, D?) + £8 (EC, t), EC, t)), 


[asf | slee Z le- Dharo NE 








< Koz | (P+ IoP +K [ KES, 


t 
f |5| < KllVcllLw (0;L~ (0)) y || (u = U) Ollia : [VEC ¢)]| 


< K( eGD? + lo, t2) + ell VEGON’, 


t t 
fasz] lu — U(rz(0))2 - 
0 0 








0 t t 
vo] IES K f (lol? + toi?) +x [iver 
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第 3 期 
根据 嵌入 定理 和 归纳 假设 (11)， 可 得 不 等 式 


Iz Ot rle 2(J;L? ae = K|% 3 加 2(J;L°(Q)) 一 


于 是 


t 
9 
[wis a (| - Dlana tle- Dlan Ulee 


AZ z| 2(J; O IN 2 C3 cz 


<i EO) + KP +o.) + f 2P, 
AES (AC Re | lemon 
OU 
jz 
<x f vel + K (lo, OP + lot DI? rf Eak 


因为 E(xz,0) = 0， 故 


fades vev = 5 2 (DIU)VE- VE) (t) - [ = ADU) oe. ve. 


I|Vell z(4;22(9)) ° IV Ellz2(7;2(0)) 








L2(J; aes 


BAO U)VE - VE) (t) > (dm + olU)|) ) D fve] 


EETh EET 


<K ‘i IVEI. 











P 


Hoy, BIEK, IRAK: a or。 联 立 上 述 各 估计 式 ， 便 得 浓度 方程 的 误差 估计 
Í “PSI + vec DIP + EGOI + JEEG EGD) 
< K I (iol? + el + toi? + uci? + neie + vel? + ver + IE R) 
+x|ve OP + Kol? + Kalot OF + | KED 
t 
af g B(iveme |e al) ee +) 


+K D (Że a alb + Heal )- (17) 
yern 7 
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接 下 来 导出 耦合 问题 的 误差 估计 . 
定理 1 SBR), po DWAR Re, 2 ApHEMS, Az, jE 和 wp 分别 为 它们 在 单 
TETKE. 整数 rp 表示 单元 马上 间断 有 限 元 空间 的 阶 ，k 表 示 RT, TARN. 假定 





min(rg, àg — l, yg — 1,wg —1,k+1)>1, VEET, (18) 


设 (pu, c) 为 问题 (1) 的 解 ， 满 足 


pEL (Hh), CEPET), So e (IH). 


假设 p, Vp, cH Vc AURA RH, id (P,U,C) 为 耦合 格式 (2)-(5) 的 解 . ABA 


| O(c — C) 
ot 











LAUJ?) + lle — Cllze( ma) 


+(J§ ((e- C(t), (e — O) 4)))? + |u- 0) 


< 天 >》， 一 








ECT, 
min(rg,Ag—1) hamle, LE—1) 
+K SO (Ea lest 2a 
pre 3/2 E> phe 3/2 fog, E i 
已 ET E E 


证 明 注意 到 对 任意 的 w v(0) = 0， 有 


lor = f Shoo? = [alot - [2] <。 [le al xf lel. 


由 于 &(:,0) =0, EER PR o = 上 <， 将 式 (14) HU Ko + 1， 关 于 时 间 积 分 后 加 到 式 (17) 中 ， 
并 利用 文献 [6] 中 的 迹 不 等 式 引 理 3.1 和 逆 不 等 式 引 理 3.2、jp 通 近 性 质 、 插 值 投影 误差 估计 
和 Gronwall 不 等 式 ， 得 到 


le D+ [ Eh) + vec ol + ecol + EGNE” 


ine wpE—1) 
< K D kee- ~ kwes—l/2 lpllws, E 


五 ET 
ad 入 已 一 1) le LE—1) 
HÈ (a rae 3/2 lelen t rhe 3/2 | me (19) 
ETh 








SAH p, n 和 (5 的 估计 式 以 及 三 角 不 等 式 ， 即 可 证 得 结论 . 


5 “归纳 假设 的 证 明 
接 下 来 检验 归纳 假设 (10)-(11). 根据 假设 (18) 和 误差 估计 (19)， 并 运用 逆 不 等 式 lu|。.o < 
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Kho |u|, HERH eV, WH 


llolIz=(4;22°(a) < Kholer) 


min(k+1,wg—1) 


= h 
< Kh! 》， ent Plone 
ECT, 








i hese Ane E,UE—1) 
+Kh( D (Exell s+ 2 ay) 


TURE h EMEB olreum) < L BRIER T zt (10). 为 证 明 归纳 假设 (11)， 首 先 
对 式 (12) 和 式 (13) 关于 时 间 t 求 导 并 取 w =r, v= 2, 48 
Oo Oo O((a(C) — a(c))u) Oo ‘a(C ) 3C Oa 
CO a? a 27) - (ae a” Hy) tY 
ða(C) ƏC _ alc) Ac Ou do 
a= |(a aa aC BB ðe me (a(C) — a(c)) a ai) 





cf Gel + «(ISL +L ter +r) 
运用 归纳 假设 (10)， 可 以 得 出 
nl (AOE - Be 





OC \dt at Ot 
< KF |(ot + letwa: |5) <el gel +e + | Sel 
lael <x (lal + lael +e + ter? + eie) +2 


注意 到 c(.,0) = 0, H Gronwall 引 理 和 估计 式 (19)， 推 出 


Ke Fl. 
IZ lee 2(J; may < L2(L?(Q)) 


k+l,wg— 1) 
h™ min( 
E 


< Ka 》， orai llles E 
EETh 


min(7r 已 ,入 已 一 in(rg,uE—1) 


+Kh 1 (Salone BF 
h 








a): 


当 九 充分 小 时 ， 可 以 得 到 


<L. 
IF Ot ae 2(J;L°°(Q)) 
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Error Analysis of a Combined Mixed Finite Element and 
Discontinuous Galerkin Method for Incompressible Miscible 
Displacement Problems 
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(College of Science, Hunan Institute of Engineering, Xiangtan 411104) 


Abstract: The numerical modeling of the miscible displacement in porous media is important and 


interesting in the oil recovery field, the environmental pollution problem and so on. A combined ap- 


proximation is applied in this paper to the given miscible displacement problem including molecular 


diffusion and dispersion. The mixed finite element method is used for the pressure equation, and the 


concentration one is solved by the non-symmetric interior penalty discontinuous Galerkin method. In- 


duction hypotheses are used to avoid the inconvenience of the cut-off operator. Based on the properties 


of the interpolation projection, an a priori hp error estimate is obtained. 
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